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Abstract— Sources  of  enhanced  ionization  in  the  Space  Shuttle  environment  are  considered,  with  particular 
emphasis  on  the  plausibility  of  a  plasma  discharge  resulting  from  the  Critical  Ionization  Velocity  (CIV) 
mechanism.  These  sources  are  then  compared  with  ion  loss  processes.  It  is  found  that  within  the  shuttle 
environment  the  loss  of  ions  is  faster  than  ion  production  from  CIV,  leading  to  the  conclusion  that  the 
observed  enhancement  in  ionization  arises  from  sources  other  than  plasma  interactions.  It  is  suggested  that 
the  seemingly  anomalous  ionization  levels  can  be  attributed  to  firings  of  thrusters  of  the  attitude  control 
system.  In  addition,  the  connection  between  the  shuttle  glow  phenomenon  and  the  enhanced  ionization  is 
considered  and  found  to  be  circumstantial  and  not  causal,  thereby  strengthening  the  argument  that  the 
primary  source  of  the  shuttle  glow  is  chemical  interactions  of  atmospheric  gases  with  surfaces  or  with 
surface-adsorbed  gases. 


INTRODUCTION 

Measurements  of  the  environment  near  the  Space 
Shuttle  have  indicated  the  presence  of  two  unexpected 
phenomena:  the  shuttle  glow  (Banks  el  al .,  1983; 
Mende  el  at.,  1983)  and  regions  of  enhanced  ion¬ 
ization  (Shawhan  and  Murphy,  1983;  McMahon  el 
al.,  1983;  Siskind  el  al.,  1984;  Raitt  el  at.,  1984; 
Shawhan  el  al.,  1984).  A  number  of  theories  have 
been  proposed  to  explain  the  origin  of  the  shuttle 
glow.  These  theories  and  the  shuttle  glow  phenom¬ 
enon  have  been  reviewed  in  several  recent  publications 
(Green  et  at.,  1985;  Mende  and  Swenson,  1985;  Kof- 
sky  and  Barrett,  1986;  Murad,  1987a,  b).  All  but 
one  of  the  theories  have  invoked  traditional  chemical 
reactions,  either  gas  phase  or  gas-surface  reactions, 
which  did  not  depend  on  the  charged  particle  back¬ 
ground.  The  one  exception  has  been  the  plasma  theory 
of  Papadopoulos  (1984),  which  has  made  an  associ¬ 
ation  between  the  shuttle  glow  phenomenon  and 
reports  of  enhanced  ionization  in  the  shuttle  environ¬ 
ment.  In  making  that  association,  Papadopoulos 
(1984)  invoked  the  Critical  Ionization  Velocity  (CIV) 
mechanism  as  a  source  of  the  enhanced  ionization,  a 
consequence  of  which  would  be  the  excitation  of  ions 
to  upper  electronic  states.  These  excited  molecular 
and  atomic  ions  would  subsequently  radiatively  de- 
excite  in  the  atmosphere  about  the  shuttle.  On  the 
other  hand,  the  most  likely  chemical  explanation  for 
the  origin  of  the  glow  is  the  excitation  of  NOj  by 
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the  surface-catalyzed  recombination  of  O  and  NO  on  K  “j 
shuttle  surfaces  (Swenson  el  al.,  1985).  Supporting 
this  idea  are  the  mass  spectrometer  data  from  theri^J 
Dynamic  Explorer  satellite  (Engebretsen  and  Hedin, 

1986)  and  from  the  Space  Shuttle  (von  Zahn  and  , 
Murad.  1986)  which  suggest  that,  indeed.  NO,  is  pre-  * 
sent  in  sufficient  quantities  to  be  consistent  with  this 
mechanism.  * 

We  explore  in  this  report  the  causal  connection 
between  the  observations  of  enhanced  ionization  and 
the  shuttle  glow  and  examine  in  more  detail  the  appli-  g 
cability  of  the  CIV  hypothesis  as  a  source  of  enhanced 
ionization  in  the  shuttle  environment.  Although 
others  have  questioned  the  connection  from  an  obser¬ 
vational  standpoint  (Kofsky,  1984;  Yee  et  al..  1984; 
Slanger,  1986),  we  consider  here  the  fundamental  pro-  ' 
cesses  which  sustain  and  drive  CIV  discharges,  and  ,  a 
examine  the  likelihood  of  its  occurrence  on  the  shuttle.  * 

The  CIV  theory,  which  was  used  by  Papadopoulos, 
was  first  suggested  by  Alfven(  1954)  as  the  cornerstone  ; 
of  his,  theory  for  the  formation  of  the  solar  system. 

The  basis  of  his  argument  was  that  a  neutral  gas, 
travelling  transverse  to  a  magnetized  plasma  at  a  rela¬ 
tive  velocity  a  critical  velocity  r*.  is  ionized,  r*  is 
given  by  \2.9^(IPIM),  where  13.9  is  a  conversion 
constant,  IP  is  the  ionization  potential  of  the  neutral 
atom  in  electronvolts,  and  M  the  mass  in  amu.  A  series 
of  laboratory  experiments  using  homopolar  devices 
(Alfven,  I960;  Fahleson,  1961)  and  neutral  beams 
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(Danielsson,  1970.  1973)  have  provided  strong  evi¬ 
dence  for  the  existence  of  CIV.  Reviews  of  the  experi¬ 
mental  evidence  and  of  the  theory  have  recently  been 
published  by  Danielsson  (1973)  and  Newell  (1985), 
respectively. 

This  report  is  divided  into  several  sections  which 
detail  the  requirements  for  the  discharge  to  proceed 
within  the  finite  geometry  of  the  shuttle  and  in  the 
presence  of  plasma  loss  processes.  An  alternate  mech¬ 
anism  which  could  explain  the  observations  of  dis¬ 
turbed  plasma  regions  near  the  shuttle  in  the  absence 
of  CIV  is  proposed.  Finally,  an  estimate  of  the  glow 
intensity  which  is  based  solely  on  reported  measure¬ 
ments  of  the  disturbed  regions  is  compared  to  the 
measured  intensity  of  the  glow  and  discussed. 

DISCUSSION 

(A)  Plasma  considerations 

(i)  Electron  impact  ionization  times.  It  is  implicit  in 
the  CIV  theory  that  ionization  is  initiated  and  propa¬ 
gated  by  electron  impact.  Since  the  energy  distribution 
of  electrons  in  the  ionosphere  is  thermal  in  nature, 
their  energies  are  below  the  ionization  potential  of 
all  the  elements.  Thus,  it  is  necessary  to  generate 
high  energy  tails  if  enhanced  ionization  is  to  occur 
in  the  vicinity  of  the  shuttle.  The  results  reported  by 
McMahon  et  at.  (1983)  indicate  that  there  is  an  en¬ 
hancement  in  electron  densities  for  electrons  in  the 
range  20  1 00  eV.  We  consider  here  the  case  of  an  elec¬ 
tron  having  an  average  energy  of  50  eV,  at  which 
value  the  electron  impact  ionization  cross-section  is 
about  I  ( —  1 7)t  cm2  for  gaseous  species  such  as  N2. 
although  for  metals  it  can  be  appreciably  higher 
(Mark,  1985).  The  ionization  time,  Tj,  is  given  by : 

t,  =  l/[/V-<x(£)-F(£)]  (I) 

where  iV  is  the  neutral  density,  <r(E)  the  electron 
impact  ionization  cross-section  as  a  function  of 
energy,  and  V  the  velocity  of  a  hot  electron.  For  a 
neutral  concentration  of  ~  1(10)  cm  '  [cor¬ 
responding  to  a  pressure  of  3(  —  7)  torr],  a(E)  of 
1  ( — 17)  cm2,  and  a  velocity  of  1(8)  cm  s  ',  the  ion¬ 
ization  time  turns  out  to  be  about  0.1  s,  thereby 
placing  an  upper  limit  on  the  time  scale  of  the  loss 
mechanisms.  If  the  losses  occur  at  shorter  times, 
then  the  discharge  can  not  be  sustained. 

(ii)  Townsend’s  criterion.  In  order  for  the  CIV 
mechanism  to  proceed,  there  must  exist  a  “seed" 


t  Throughout  this  paper,  a(b)  denotes  a  x  10*. 


density  of  energetic  electrons  within  the  discharge 
region  capable  of  ionizing  the  background  gas.  The 
Townsend  criterion  stipulates  that  in  order  to  sustain 
a  discharge,  there  must  be  on  average  at  least  one 
electron  impact  ionization  during  the  lifetime,  x,  of  an 
electron  within  the  discharge  region.  This  implies 
that  an  electron  heated  by  a  plasma  process  (as  yet 
undefined)  has  to  have  enough  time  to  ionize  another 
neutral  atom  or  molecule  during  its  transit. 

Within  the  theory  of  quasilincar  plasma  physics 
(Formisano  et  al..  1982),  it  is  calculated  that  thermal 
electrons  can  absorb  electrostatic  wave  energy  to  form 
a  high  energy  tail  or  plateau  in  their  energy  distri¬ 
bution.  The  formation  of  this  hot  tail,  rather  than 
bulk  electron  heating,  renders  the  heated  electrons 
energetic  enough  to  cause  ionization  of  the  neutral 
gas.  Only  a  small  fraction,  *,  of  the  total  population 
of  electrons  migrates  into  the  tail,  while  the  bulk  of 
the  electron  population  is  essentially  unhealed,  as 
illustrated  in  Fig.  I.  The  unhealed  electrons,  in  con¬ 
trast  to  those  in  the  tail,  have  insufficient  energy  to 
ionize  the  neutral  gas.  Consequently,  the  Townsend 
criterion  needs  to  be  modified  so  that  a  heated  electron 
has  to  make  at  least  a  1  ionizations  during  its  lifetime. 
X/  .  within  the  discharge  region.  Mathematically,  the 
modified  Townsend  criterion  becomes : 

a  '[t,+x„]<x,  (2) 

where  x,  is  the  ionization  lime,  and  L  is  the  transverse 
dimension  of  the  discharge  region.  The  lime.  x,.  of  a 
heated  electron  in  the  discharge  region  is  given  by : 

x,  =  L/V.  (3) 

Tanaka  and  Papadopoulos  (1982)  obtained  the  fol¬ 
lowing  expression  for  the  electron  heating  time : 

rH  -  30/m,  „  (4) 

where  to,  „  is  the  lower  hybrid  frequency. 

In  the  shuttle  environment,  equation  (4)  typically 
yields  x„  -  1  10  ms,  depending  on  the  ion  mass  and 
the  local  magnetic  field  strength.  To  calculate  the 
inequality  in  equation  (2),  we  use  the  following  typical 


Fig.  I.  Electron  energy  DtsTRn»i;xxoN  resulting  from 
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parameters:  <£>  %  5-30  eV  (Papadopoulos,  1984), 
L  ~  6  m,  and  a  ~  0.01,  JV  ~  1  (10)  cm  ',  and 
<r~0.5(— 16)  cnr.  As  a  result  inequality  (2)  gives 
about  2  s  on  the  left  side  and  about  6  /ts  on  the  right : 
the  two  sides  differ  by  orders  of  magnitude.  Thus 
Townsend’s  criterion  renders  a  sustained  discharge 
unlikely  in  this  environment. 

In  order  to  examine  the  parametric  dependence  of 
the  modified  Townsend  criterion,  equation  (2),  it  is 
convenient  to  define  a  ratio  F: 

F(N.V,L)  =  x  '( t,+th)/t,,  (5) 

When  f(jV)  <  I,  equation  (2)  would  be  satisfied.  Fig¬ 
ure  2  shows  the  dependence  of  F(N)  on  the  neutral 
density,  N ,  for  various  values  of  V  starting  from 
that  of  ion  sound  velocity  [F  ~  Cs  ~  1(5)  cm  s  1 : 
Papadopoulos.  1984],  while  the  parameter  L  is  taken 
as  6  m.  F(N)  is  a  monotonically-decreasing  function, 
because  the  probability  of  ionization  increases  as  N 
increases.  The  slope  of  F(N)  decreases  to  zero 
asymptotically,  because  at  high  N  t,  becomes  short 
compared  to  the  heating  time,  tH,  so  that  the  latter 
gradually  becomes  dominant.  As  a  function  of  spatial 
dimension,  L,  In  F(L)  decreases  linearly  with  In  L 
(Fig.  3).  In  this  figure,  F(L)  is  plotted  with  N  ~  l(!0)cm  ' 
(Narcisi  et  al.,  1983)  and  the  same  value  of  V  as  in 
Fig.  2.  F(L)  is  greater  than  unity  until  L  becomes 
quite  large.  This  implies  that  CIV  discharge  can  be  sus¬ 
tained  only  when  L  is  very  large  (perhaps  -  I  km). 
The  shuttle  tail  is  only  6  m.  Since  the  energy  range 
of  a  hot  electron  tail  in  the  shuttle  vicinity  is  esti¬ 
mated  to  be  5  30  eV  (Papadopoulos,  1984)  the  velo¬ 
city  V  of  a  hot  electron  traversing  a  discharge  region 
is  expected  to  be  much  greater  than  Cs.  With  these 
parametric  values  we  conclude  that  F  »  1 .  A  higher 
neutral  density,  for  example  I  ( 1 1 )  cm  \  would  not 
change  the  conclusion.  Therefore  a  discharge  in  the 


Fig.  2.  Variation  of  F  (see  text  for  definition)  as  a 

FUNCTION  OF  NEUTRAL  NUMBER  DENSITY. 


Fit;.  3.  Variation  of  F  as  a  function  of  confinement 
LENGTH. 

See  text. 


shuttle  environment  is  highly  unlikely. 


(iii)  The  wave  vector  triangle.  The  relevance  of  the 
lower  hybrid  (LH)  wave  heating  of  electrons  to  the 
CIV  process  has  been  studied  by  several  authors 
(Raadu.  l978;Galecv,  1981 :  Formisano  et  al..  1982: 
Papadopoulos,  1984,  1985;  Abe  and  Machida,  1985; 
Machida  and  Gocrtz.  1986).  Basic  investigations  of 
the  LH  wave  electron  heating  have  been  performed 
by  McBride  et  al.  (1972).  McBride  et  al.  (1972)  and 
Papadopoulos  (1984)  have  calculated  that  the 
maximum  efficiency  of  heating  occurs  at  the  wave 
vector  angle  0  =  I .  where  0  is  given  by : 

0  =  (kJkHM/m)' :.  (6) 

The  wave  vector  k  is  related  to  the  wavelength. 
and  to  the  hybrid  frequency,  w,  H.  by : 

k  =  In!/.  (7) 

k  =  wLH/r  (8) 

where  v  is  the  beam  velocity  of  the  ions  energizing  the 
waves. 

0  =  I  leads  to  a  requirement  that  the  size  of  the 
wave  vector  triangle  (Fig.  4)  be : 

A.  1  -  (\/k)(Mlm)'2.  (9) 


The  size  requirement  for  CIV  in  space  has  been  dis¬ 
cussed  by  Newell  (1985)  and  Kelley  et  at.  (1986) 
in  connection  with  the  seemingly  different  observed 
efficiencies  for  CIV  in  the  laboratory  (high  efficiency) 
and  in  space  (low  efficiency).  We  examine  below  the 
size  requirement  (equation  (9)]  for  the  case  of  homo- 
polar  devices  and  for  the  case  of  space  experiments.  p0(jag 
From  equation  (4)  and  taking  the  electron  heating  — 
time,  tH,  to  be  3  /is  (Fahleson,  1962)  we  obtain :  i/or 
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Fig.  4.  Wave  vector  triangle. 


<0i.h  ~  1(4)  kHz.  (10) 

Substituting  this  value  of  o»LH  into  equation  (9).  and 
assuming  v  ~  37  km  s  '  for  H,  (Fahleson,  1962),  we 
obtain  /  ~  14  cm,  which  is  smaller  than  the  dimen¬ 
sions  of  a  homopolar  device  (Fahleson,  1962,  gives 
dimensions  of  IS  cm  height  x  20  cm  radius). 

For  the  Space  Shuttle  conditions,  substituting 
«/LH  -v  3.6  kHz  (Kelley  et  a/.,  1986)  and  v  ~  2x7.3 
km  s  '.we  obtain  /  ~  1  km,  which  agrees  with  Kelley 
et  al.  (1986).  This  is  about  two  orders  of  magnitude 
larger  than  the  Space  Shuttle  tail  cross-section.  In 
other  words,  the  wave  vector  triangle  does  not  fit  into 
the  transverse  direction  of  the  ion  beam. 

(iv)  Wave  vector  condition  for  a  finite  system.  We 
consider  here  the  case  of  a  system  with  finite  dimen¬ 
sions  having  a  simple  configuration  as  shown  in  Fig. 
5a :  the  ambient  magnetic  field  B  is  parallel  to  the  z- 
axis  and  the  ion  beam  parallel  to  the  .v-direction.  The 
stability  analysis  seeks  a  solution  of  the  form 

<p  =  <p(z)exp(  —  uot+ikx)  (II) 

resulting  in  an  eigenvalue  problem  of  the  form  (see 
Appendix) 

(d/dr)  {*.  (r)(d^o(c)/dr) }  +  T2<p(z)  =  0  (12) 

where 

T2  =  —k2r.fk,  at),  (13) 

<p  is  the  plasma  wave  potential,  k  the  wave  vector 
( =  kx  in  this  geometry),  to  the  wave  frequency,  and  r, 
the  dielectric  function. 

The  eigenvalue  equation,  (13),  can  be  solved  for 
given  simple  density  profiles.  We  consider  a  square 
density  profile  (Fig.  5b),  in  which  the  plasma  densities 
are  n  inside  the  beam  (Region  I)  and  %  outside  the 
beam  (Regions  II  and  III),  while  the  boundaries  of 


Fig.  Sa.  Coordinates  for  a  simple  model  of  a  finite-sized 

BEAM  PROPAGATING  ALONG  THE  ^-DIRECTION  PERPENDICULAR 
TO  THE  AMBIENT  MAGNETIC  FIELD  B.  WHICH  IS  IN  THE  r- 
DIRECTION. 
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Fig.  5b.  A  simple  density  profile. 


the  beam  are  r  =  ±L.  For  this  case  equation  (12) 
reduces  to  (see  Appendix) 

dV/dr-r>  =  0  (14) 

where 

r2{k,u>)  =  k2r.x(k,ti))  lc.(k,(o).  (15) 

The  appropriate  boundary  conditions  are  <p  being 
finite  at  infinity  and  ip  as  well  as  d<p/dz  being  con¬ 
tinuous  at  the  beam  boundaries.  Aside  from  the  trivial 
solution  (p  =  0.  we  obtain,  as  shown  in  the  Appendix, 
the  condition 


T  tanh  (TZ,)  =  — T0  (16) 

where  T0  is  the  function  T  in  Region  II  (or  III),  where 
the  density  is  «,>•  Using  this  simple  density  profile  we 
have  solved  equation  (16)  numerically:  the  model 
assumes  a  finite  kt  and  no  k,  [see  equations  (Al) 
and  (A2)  in  the  Appendix],  Numerical  results  are 
shown  in  Figs.  6a  and  6b.  In  the  infinite  beam  limit 
(kL  -» oo),  the  maximum  growth  rate  (Im  «>/«>,  - 
0.19)  occurs  at  kr /w,  -  0.8  (Fig.  6a).  Figure  6b 
shows  the  finite  beam  case  where  the  growth  rate 
(Imru/tu,)  at  kv/cji  ~  0.8  is  plotted  as  a  function 
of  kL ;  Im  w/a>|  reaches  an  asymptotic  value  of  0. 19  as 
kL  — f  cc,  which  agrees  with  the  infinite  beam  result  in 
Fig.  6a.  It  is  interesting  to  observe  that  in  Fig.  6b 
the  instability  growth  rate.  Imw/iu„  decreases  with 
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Fig.  6a.  Dispersion  relation  of  a  simple  model  of  an 

INFINITE  BEAM  PROPAGATING  PERPENDICULAR  TO  THE  AMBIENT 
MAGNETIC  FIELD. 

The  imaginary  root  lm<o  corresponds  to  the  Re  to  branch 
that  starts  at  to  =  0.  Max  |Ima>  «j,]  ~  0.19  occurs  at  about 
ki  lto,  ~  0.8.  The  parameters  used  are  fl Jto,  =  I.  and 
njn„  =  0.1.  The  ions  are  O.' . 


Fkl  6b.  Numerical  solutions  of  Equation  (16)  as 

FUNCTIONS  OF  FINITE  kL. 

The  parameters  used  are  the  same  as  in  Fig.  6a,  with 
ki  /to,  =  0.8.  Notice  that  the  growth  rate,  1m  <o/<o„  decreases 
as  Art  decreases,  and  that  in  the  Art  -»  oo  limit,  one  recovers 
the  maximum  growth  rate  Imw/ru,  ~  0.19.  agreeing  with  the 
results  in  Fig.  6a. 


decreasing  kL.  For  example,  taking  a  shuttle  tail 
dimension  of  2L  -*•  6  m  and  a  wavelength,  i,  '6  m 
(i.e.  kL  ^  it),  the  instability  growth  rate  is  reduced  to 


~  1/5  of  that  in  the  infinite  plasma  case.  For  longer 
wavelengths,  the  growth  rate  would  be  even  lower. 
This  result  suggests  that  long  wavelength  plasma 
instabilities  obtained  using  infinite  plasma  theories 
may  be  rendered  inappropriate  (i.e.  inefficient)  by  the 
finite  spatial  extent  of  the  surfaces. 

(B)  Ion  loss  mechanism 

(i)  Diffusion.  Ionization  generated  by  the  CIV  pro¬ 
cess  in  the  vicinity  of  the  Space  Shuttle  can  be  lost  by 
diffusion.  If  the  diffusion  time  is  shorter  than  the 
ionization  time,  then  the  plasma  density  cannot  be 
increased.  Parallel  to  a  magnetic  field,  diffusion  of  a 
plasma  with  a  neutral  gas  background  is  governed 
by  the  ambipolar  diffusion  process.  The  ambipolar 
diffusion  coefficient,  Da,  is  ~  2D,,  where  D,  is  given 
by: 

D,  =  k T,/Mv  =  k T'  [M •  N ■  om  •  r]  (17) 

where  M.  N.  r.  and  <rm  arc  the  ion  mass,  number 
density,  velocity,  and  ion  neutral  collision  cross- 
section.  respectively.  Across  a  magnetic  field.  B. 
diffusion  of  a  plasma  is  given,  approximately,  by 
Bohm's  semi-empirical  formula  (Chen,  1974): 

Dh  =  krc/(!6-p-B).  (18) 

In  order  to  estimate  the  magnitude  of  this  loss  pro¬ 
cess,  we  consider  the  case  where  k  Tr  %  30  eV.  B  s;  0.3 
G,  M  ~  16,  and  o,„  %  I ( —  1 5)  enr  (Hamilton  and 
Knewstubb,  1984).  The  scale  length  is  taken  as  6  m. 
which  is  approximately  the  length  of  the  shuttle  tail 
as  well  as  the  gyroradius  of  an  ion  with  about  the 
shuttle  speed.  For  a  simple  cylindrical  geometry,  the 
diffusion  time,  rd  =  R:/{2D).  where  (Chen.  1974) : 

D=D.  +  DH.  (19) 

Figure  7  shows  the  competition  between  diffusion  time 
and  electron  impact  ionization  time.  If  we  equate  r,  to 
rd,  we  obtain  a  critical  value  for  N.  Nc  (the  crossing 
point  in  Fig.  7),  above  which  ionizition  dominates 
over  diffusion,  and  below  which  diffusion  dominates. 
This  value  of  JVC  turns  out  to  be  -  5(9)  cm  \ 

(ii)  Dissociative  electron  ion  recombination.  Once 
CIV  has  occurred,  the  discharge  will  be  quenched  if 
the  electron  and  ion  densities  exceed  a  value  whose 
magnitude  depends  on  the  recombination  coefficient. 
Radiative  recombination  of  atomic  ions : 

A*  +e-»  A+hv  (20) 

where  A  is  an  atomic  ion.  is  exceedingly  slow.  A*,  being 
~4(  — 12)  cm’  molec  1  s  1  (Bates  and  Dalgamo. 
1962 ;  Dunn  etal.,  1984).  However,  the  rate  coefficient 
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Fig.  7.  Competition  between  ionization  and  diffusion 
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for  dissociative  recombination : 

AB*  +  e- A  +  B  (21) 

where  AB  is  a  generic  molecular  ion,  is  much  larger. 
k2,  being  of  the  order  of  l(—  6)  to  !(  —  7)  cm'  molcc  1 
s  '  (Bates.  1979;  McGowan  and  Mitchell,  1984).  The 
increased  pressure  near  the  shuttle  is  due  to  molecular 
contaminants  such  as  H.O  and  CO:  (Wulf  and  von 
Zahn,  1986),  and  measurements  suggest  that  ions 
from  these  molecules  (Narcisi  el  al.,  1 983 ;  Grebowsky 
el  al.,  1988)  are  the  primary  sources  of  enhanced 
ionization.  Loss  by  electron  ion  recombination  would 
involve  these  contaminants,  as  discussed  recently  (Murad 
and  Lai,  1986).  Briefly,  there  are  three  immediate 
consequences  of  this  loss  mechanism:  (I)  the  density 
of  electrons,  crucial  in  sustaining  the  chain  ionization 
process,  is  reduced  by  recombination;  (2)  consider¬ 
able  energy  loss  by  hot  electrons  is  lost  as  a  result  of 
radiation  liberated  in  reaction  (21);  and  (3)  a  lower¬ 
ing  of  the  mass  of  the  neutral  species  would  result 
in  a  corresponding  increase  of  the  critical  velocity 
threshold  discussed  earlier  (Axnas,  1978).  Depending 
on  the  parametric  conditions,  CIV  may  proceed  in  an 
oscillatory  manner,  or  it  may  be  completely  quenched 
(Lai  el  al.,  work  in  progress).  If  the  relative  velocity 
between  the  neutral  gas  and  the  plasma  is  not  high 
enough  to  replenish  the  losses,  the  discharge  process 
would  be  quenched.  In  the  case  of  the  Space  Shuttle, 
the  energy  budget  of  CIV  is  already  marginal ;  it  relies 
on  the  8%  ambient  ions  (Papadopoulos,  1983,  1984) 
reflected  on  the  ram  side  of  the  Space  Shuttle.  Any 
loss  mechanism  would  make  the  initiation  of  the  CIV 
discharge  even  more  difficult. 


(iii)  Line  excitation.  Neutral  species  can  be  excited 
to  allowed  electronic  states  by  electron  bombardment, 
which  would,  in  turn,  radiate  very  fast  (within  a  few 
nanoseconds).  The  net  effect  of  such  a  process  is  to 
cool  the  electrons;  Newell  and  Torbert  (1985)  have 
calculated  that  nearly  50%  of  the  electron  energy  is 
lost  in  neutral  heating  in  the  shuttle  environment.  This 
becomes  a  substantial  loss  mechanism  which  would 
adversely  affect  an  already-marginal  CIV  energy 
budget,  as  discussed  above.  A  loss  of  50%  of  the 
energy  would  tend  to  quench  the  nascent  discharge. 
A  more  detailed  study  of  neutral  excitation  and  de¬ 
excitation  processes  in  molecular  CIV  systems  is  in 
progress  (Murad.  1987a ;  Lai  el  al.,  1988). 

CONFUSIONS 

(A)  Ionization 

The  origin  of  the  observed  enhancement  in  ion¬ 
ization  in  the  vicinity  of  the  Space  Shuttle  is  prob¬ 
lematic.  since,  as  discussed  above.  CIV  is  unlikely  to 
be  sustained  in  the  vicinity  of  the  Space  Shuttle.  This 
conclusion  then  requires  an  alternative  explanation 
for  the  origin  of  ionization  in  the  vicinity  of  the  Space 
Shuttle.  A  likely  source  is  from  thruster  firings,  since 
mass  spectrometric  measurements  in  the  bay  of  the 
Space  Shuttle  indicate  a  direct  correlation  between 
densities  of  ions  (Narcisi  el  al..  1983:  Grebowsky  el 
al..  1987)  and  neutrals  (Narcisi  el  a/.,  1983;  Wulf  and 
von  Zahn.  1986)  with  thruster  firings.  While  the  flow 
fields  from  thrusters  to  the  shuttle  bay  are  not  known 
at  this  time,  it  seems  reasonable  to  attribute  the 
observed  enhancement  in  ionization  level  to  large 
amounts  of  ions  released  in  the  exhaust  of  thruster 
engines. 

(B)  Glow 

Mcnde  and  Swenson  (1985)  have  estimated  the 
glow  intensity  to  be  10  100  kR  [I  R  =  1(6)  photons 
cm  '  s'].  This  corresponds  to  a  photon  emission  rale 
of  I  (10)- 1(1 1)  photons  cm  ~s  '  per  column  length, 
this  being  the  length  of  the  tail  section  of  the  Space 
Shuttle  (about  6  m).  In  the  plasma  mechanism  the 
ionization  involves  electron  impact  ionization  and 
excitation  to  upper  electronic  states  of  the  ions : 

cf+AB^{AB*)*  +  2e  (22) 

[AB'  )*  -+AB*  +Av.  (23) 

AB  is  any  neutral  species  (atomic  or  molecular),  e,  a 
fast  electron,  and  {AB*  )*  an  excited  state  of  AB' . 
The  rate  of  formation  of  {AB*)*  is  given  by 
ff 2 2 *  ^'[crl'M/l),  where  a21  is  the  cross-section  for 
reaction  22;  V  the  velocity;  [cf),  the  concentration 
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of  fast  electrons;  and  [AB]  the  concentration  of  the 
neutral  species  AB.  For  CO,  at  an  electron  energy  of 
1 50  eV,  <t,_,  is  -  0.5(  —  16)  cm:  (Mental!  el  al.,  1973). 
At  20-30  cV.  ffv,  is  likely  to  be  much  smaller;  a  con¬ 
servative  estimate  would  be  1  ( —  17)  cm2.  For  [er],  we 
use  a  value  of  2(5)  cm  \  as  has  been  reported  by 
McMahon  el  al.  ( 1983),  and  for  [AB]  we  use  a  value 
of  1(10)  cm  ’.  Using  1(8)  cm  s  we  calculate 
the  rate  of  formation  of  (AB*)*  to  be  2(7)  cm  's  '. 
For  a  path  length  of  6  m.  this  density  corresponds  to 
a  column  density  of  1 ,2(  10)  molecules  cm  2  s  *.  If  the 
emission  rate  is  high,  e.g.  10%,  then  the  predicted 
intensity  would  have  a  maximum  of  1.2(9)  photons 
cm  2  s  ',  or  1  kR.  This  corresponds  to  0.1-0.01  of 
the  observed  intensity.  Thus  we  conclude  that  the 
plasma  hypothesis  cannot  explain  the  shuttle  glow 
phenomenon. 

Yee  et  al.  (1984)  reached  a  similar  conclusion  based 
on  the  lack  of  correlation  of  shuttle  glow  intensity 
with  plasma  density.  More  recently,  Slanger  (1986) 
found  a  lack  of  correlation  with  neutral  density  in  the 
V  section  of  the  tail.  Moreover,  the  near-continuum 
nature  of  the  glow  (Swenson  et  al..  1985)  has  lent 
credence  to  the  suggestion  that  the  glow  is  caused  by 
NO,  formed  by  recombination  on  the  surfaces  of  the 
shuttle.  Indeed,  mass  spcctromctric  data  support  this 
finding  (von  Zahn  and  Murad.  1986).  The  nature  of 
the  spectrum  also  argues  against  the  plasma  hypothe¬ 
sis.  as  discussed  in  some  detail  by  Kofsky  (1984).  For 
example,  the  major  contaminant  on  the  Space  Shuttle 
is  H_,0.  The  spectrum  of  its  ion.  which  would  be 
formed  in  reaction  (22),  is  a  well  defined  band 
system  (Lew.  1976)  contrary  to  the  observation. 

(C)  Limits  on  ionization 

Figure  7  shows  a  crossing  between  diffusion  loss  of 
ions  and  generation  of  ions  by  electron  impact  at  a 
neutral  density  of  -  5  x  I0g  cm  \  This  corresponds 
to  the  density  at  an  altitude  of  -  240  km  (U.S.  Stan¬ 
dard  Atmosphere.  1976).  At  altitudes  greater  than  this 
ions  will  be  lost  by  diffusion,  while  at  lower  altitudes, 
it  is  likely  that  ionization  will  dominate  diffusion. 
At  lower  altitudes  other  loss  mechanisms,  such  as 
electron-ion  recombination,  might  become  dominant. 
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APPENDIX 

Linear  dispersion  relation 

Consider  an  infinite  slab  beam  with  density  n  and  thickness 
2 L  interacting  with  an  infinite  plasma  density  n„  in  the  pres¬ 
ence  of  a  magnetic  field  B  (Q  =  cB  mc )  in  the  r  direction, 
which  is  perpendicular  to  the  beam  propagating  with  velocity 
r  in  the  x  direction  (sec  Fig.  5a).  For  electrostatic  per¬ 
turbations  of  the  form 

i  ,=(  ,(;:) exp (-i<ur  +  ikx)  (Al) 

and 

<p  -  <p(:)  exp  ( -  wit  +  ik.v)  <A2) 

the  equation  of  motion 

Bi(dv„dr)  =  e[E,  +  ( I /'<•)(*,  x  B)J.  (A3) 

where  E,  =  —Nip.  and  the  continuity  equation 

r'n,/r/  +  V-(n0vl +n,v)  =  0  (A4) 

yield 

4nn,  =  v£k2<i>l[{(o—kr)2  —  n:] 

-  (rWr)l«»;(rty/dr))/|(»  -  A  t-]2  ( A5) 

where  the  plasma  frequency  otc  ( =  ^  (4 ntfn/m)\  is  a  function 
of  the  density,  n(z).  For  sci  I -consistency,  the  quantities  n, 
and  tp  have  to  satisfy  Gauss's  equation : 

\'tp=-4nni.  (A6) 

Eliminating  n,  from  (A5)  and  (A6).  one  obtains 

(d/dz)[e.(z)(d/dz)J<f»(z)  +  T‘>  =  0  ( A7) 

where 

c.(r)  =  I  -«>;!(<» -kr):. 
f.,(z)  =  I  — c»,?/I(fo-kr)J— QJ).  and 
r  =  -k:c,(z).  (A8) 
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To  recover  the  infinite  limit  for 

~  v>(0) exp (  - i*.r)  ( A9) 

one  obtains  from  IA7)  and  (AS)  the  dispersion  relation 
1  —o>2,iv)—kj)2{kJk)2 

~kxr)~  -02](kx/k) 2  =  0  (A10) 

where  k 2  -  k2.+k;. 

Simple  density  profile 

The  eigenvalue  equation  (A7)  can  be  solved  for  a  given 
density  profile.  For  the  simplest  case  of  perpendicular  beam 
propagation  in  a  background  plasma,  we  consider  a  square 
density  profile,  where  the  beam  is  propagating  in  the  .v- 
direction :  the  density  profile  along  the  e-direction  is  dis¬ 
played  in  Fig.  5b.  The  plasma  densities  are  assumed  to  be 
constant  inside  the  beam  (Region  1)  and  outside  the  beam 
(Regions  II  and  III).  The  beam  boundaries  are  at  :  =  L  and 
-I. 

For  this  case  of  a  square  density  profile,  the  plasma 
frequency  o>,  is  constant  (devde  =  0)  inside  each  region. 
Equation  (A7)  becomes 

d'(/j(-l  dr- r2tp(;)  =  0  (All) 

where 

=  A!e,.(r)  c.(r).  (A12) 

The  solution  of  equation  (All)  now  yields 

</»(;)  =ae  1  -t -he1'  Region!  (A  1 3) 

and 


<p(z)  =  a„e  '«■'+ A„e’»-'  Regions  II  and  III.  (AI4) 
Using  the  boundary  conditions : 

( 1 )  ip(  x )  is  finite,  and 

(2)  both  tp  and  dip/dz  are  continuous  at  ;  =  /..  /. 

we  obtain,  aside  from  the  trivial  solution  of  constant  poten¬ 
tial  which  is  not  of  interest,  the  following  dispersion  relation  : 

T  tanhlHL)  =  -T,,  (AI5) 

where  F.  defined  in  (A12),  is  a  function  in  the  beam  region 
1.  while  r„  refers  to  regions  II  and  III,  which  are  outside  the 
beam. 

Taking  into  account  both  electron  and  ion  contributions 
explicitly.  (A12)  becomes 


r: 

=  k  2e ,  (nc.  «h  )ir. ..  («„.  n, .  «h ) 

f  A 16) 

and 

r,; 

i  =  k2i:Jnc,n,.  0)/s.(wc.n1.0) 

(A17) 

where 

i ;l(ne,n1.nh)  =  1  - 

■o>;  (to2 -O;) 

— v>;  la i  2  — i  «<u  —  kr ) ;  —  Jlf  ] 

(AI8) 

and 

<:-(nc. «„»(,)  = 

1  —o)2,  or  —  o)2  or  -  n>l  U’>—kr)2. 

(AI9) 

0  refers  to  the  outside  of  the  beam.  u>e  is  the  background 
plasma  frequency,  to,  the  ion  plasma  frequency,  ft,,  is  the 
background  gyrofrequency.  B  the  ambient  magnetic  field. 

=  4ne2nJM)  the  ion  beam  plasma  frequency. 
nb  («„  =  n„+nh)  the  ion  beam  density.  the  ion  beam  gyro- 
frequency,  v  the  ion  beam  velocity,  and  k  the  wave  vector. 
The  dispersion  relation  (A  1 5)  is  derived  as  described 
previously. 


SI 


V7> 


Unclassified 


la.  REPORT  SECURITY  CLASSiFiCAfiON 

Unclassified 


4  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PACE 


Fb.  RESTRICTIVE  MARKINGS 


3  OlSTRlSUTlON  /  AVAILABILITY  OF  REPORT 


b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


ftd  M AMt  OP  i*fcKK)KMlNG  OHuAMlZATlON  6o  OFFICE  SYMBOL  7d  MAMt  OF  MOM  TO  RING  ORGANIZ  AT  ION 

Air  Force  Geophysics  Laboratory 


oc.  aODRESS  (Ory.  SMTP.  ana  ZIP  COM3 

Hanscom  AFB 

Massachusetts  01731-5000 


70.  AOORESS  (Cry,  Sum,  and  ZIP  COM3 


So.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
Of  appJrtabJa)  I 


He.  ADDRESS  (C/ty.  Sum.  and  ZIP  COM3 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 

62101F 


II  Title  (Include  Security  CJauiftcafionl 

The  Role  of  Plasma  Processes  in  Che  Space  Shuttle  Environment 


12  personal  authoh(S) 

Shu  Lai,  W.F.  Denig,  Edmond  Murad,  William  J.  McNeil* 


1 3a  TYPE  OF  REPORT  Il3b.  TIME  COVERED  Il4  DATE  OF  REPORT  lYpar,  Month.  0ay3 

Reprint  I  from _ to _  |  1988  September  6 


16  supplementary  notation  *Radex,  Inc. ,  Bedford,  MA  01730  -  Reprinted  from  Planet  Space  Sci 
Vol.  36,  No.  8,  pp.  841-849,  1988 


COSATl  COOES 


group  sub-group 


IS  SUBJECT  TERMS  (Cononu*  on  reverie  ,f  nttetury  and  identify  by  block  number! 

Space  Plasma*  Space  Shuttle,’  Critical  Ionization,’ 
Diffusion  ;  ^  ,;A  f  c  ,  t\ l>  4  )*&■ 


19  ABSTRACT  (Conbnu*  on  revets*  if  necessary  and  identity  Oy  block  number) 


^3  - 

Ahtratl— Source,  of  cnhaaonl  wnuaiwa  w  ibe  Space  SbuUk  caviroemcM  arc  cooauicml,  nub  particular 
ciupuaoa  on  the  ptauaRxni;  oF  a  Hanna  liuebargc  mubag  From  ibe  Cnucai  loauoUoa  Veka. it y  (CIV) 
inctliaataui.  1  beae  mum  ale  Darn  tompatol  nub  MB  luaa  pruccaaca.  It  u  luoml  (bat  nubia  (be  abulia 
cnaitouuwul  ilw  loaa  ol  wm  la  iaalcr  lhaA  iub  pfodiMM  From  CIV,  kri..g  to  da  .•,o^-i„^.m  dm  tkc 
.RnoiaoJ  .-i.loOK.iiwni  ta  toaiaalioa  ariaca  from  aourcea  otbcl  (baa  pi  ima  taur— .liona.  tl  a  akdawawo  (bai 
ilw  aaviunai)  anomaloua  Mauauua  leaeia  call  b»  aiinboual  lo  (map  oF  ihnuicra  at  ibe  ailuuOc  control 
’l*™  aJUinm,  Urn  coaocuw  bnneea  ibe  abuik  don  pbcnoaaaaoa  aau  (be  Cfltui.nl  -nmnrm  u 
ifluwkiui  an4  looaJ  10  6r  cuiamiublul  ami  au  tnail,  Iknb)  atfcagdlnMi  Ibe  mmol  (bat  (be 
primary  aourcc  o <  ibe  abuue  (Ion  la  irbrmaial  later  aauoafl  or  uiupbenc  gaaea  nub  am  Moca  oc  nub 
<orlacc-a4wrbc4  aaaca.  Y/  _  /) 


20.  DISTRIBUTION!  AVAILABILITY  OF  ABSTRACT  121.  ABSTRACT  SECURITY  CLASSIFICATION 

^  □uNCLASSlFlfeO/UNLIMlTEO  fl  SAME  AS  RPT.  □  QTIC  USERS  I  OnclSMifimd 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b  TELEPHONE  (JncluM  Are*  Code)  22(  OFFICE  SYMBOL 

Edmond  Murad  I  (617)  377-3176  ph* 


DO  FORM  1473, 94  ma:; 


43  APR  vSction  may  M  UMB  or  til  caFiauWeb 


88  9  AH  «1l|pTgg'tiOM  J  yot^. 


wi  vi  j  iv*;;  a 


Unclassified 


